An accurate and efficient tau-leaping procedure for the simulation of chemical reaction systems
نویسنده
چکیده
By explicitly representing well-stirred chemical reaction systems with independent, unit Poisson processes we develop a new adaptive tau-leaping procedure. The procedure developed is novel in that we enforce any leap condition via a post-leap check as opposed to performing a pre-leap tau selection. Further, we perform the post-leap check in such a way that the statistics of the sample paths generated will not be skewed by the rejection of a leap. By performing a post-leap check to ensure accuracy, the method developed in this paper is guaranteed to never produce negative population values. The efficiency of the method developed here is demonstrated on a model of a decaying dimer.
منابع مشابه
Efficient step size selection for the tau-leaping simulation method.
The tau-leaping method of simulating the stochastic time evolution of a well-stirred chemically reacting system uses a Poisson approximation to take time steps that leap over many reaction events. Theory implies that tau leaping should be accurate so long as no propensity function changes its value "significantly" during any time step tau. Presented here is an improved procedure for estimating ...
متن کامل$S$-Leaping: An adaptive, accelerated stochastic simulation algorithm, bridging $\tau$-leaping and $R$-leaping
We propose the S-leaping algorithm for the acceleration of Gillespie’s stochastic simulation algorithm that combines the advantages of the two main accelerated methods; the τ -leaping and R-leaping algorithms. These algorithms are known to be efficient under different conditions; the τ -leaping is efficient for non-stiff systems or systems with partial equilibrium, while the R-leaping performs ...
متن کاملAccurate implementation of leaping in space: the spatial partitioned-leaping algorithm.
There is a great need for accurate and efficient computational approaches that can account for both the discrete and stochastic nature of chemical interactions as well as spatial inhomogeneities and diffusion. This is particularly true in biology and nanoscale materials science, where the common assumptions of deterministic dynamics and well-mixed reaction volumes often break down. In this arti...
متن کاملEfficient Leaping Methods for Stochastic Chemical Systems
Well stirred chemical reaction systems which involve small numbers of molecules for some species have a stochastic behavior and can be modeled by a continuous time, discrete state Markov process. An exact method for simulating the time evolution of the system is the Stochastic Simulation Algorithm, but this method is extremely slow for realistic biological systems. We propose an adaptive leapin...
متن کاملAvoiding negative populations in explicit Poisson tau-leaping.
The explicit tau-leaping procedure attempts to speed up the stochastic simulation of a chemically reacting system by approximating the number of firings of each reaction channel during a chosen time increment tau as a Poisson random variable. Since the Poisson random variable can have arbitrarily large sample values, there is always the possibility that this procedure will cause one or more rea...
متن کامل